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Abstract. Compelling evidence has recently emerged from lattice simulations in favor of the massive solution of the 
Schwinger-Dyson equations of Landau-gauge QCD. The main objections to these lattice results are based on possible Gribov- 
copy effects. We recently installed at IFSC-USP a new GPU cluster dedicated to the study of Green's functions. We present 
here our point of view on the Saga and the status of our project. We also show data for the 2D case on a 2560^ lattice. 
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THE SAGA 

Intuitively, the explanation of color confinement should 
be encoded in the infrared (IR) behavior of QCD Green's 
functions. The Landau-gauge Gribov-Zwanziger (GZ) 
confinement scenario and the scaling solution obtained 
by solving Schwinger-Dyson equations (SDEs) demand 
— for any space-time dimensions D > 2 — a null gluon 
propagator at zero momentum and an IR-enhanced ghost 
propagator [1]. At present, there is wide agreement [2] 
that numerical simulations in minimal Landau gauge 
show (in the infinite-volume limit): 1) an IR-finite gluon 
propagator D(p) in D = 3,4 and a null D{0) in 2D, 2) 
violation of reflection positivity for the gluon propagator 
in D — 2,3,4 and 3) an essentially free ghost propagator 
G{p) in D = 3,4 but IR-enhanced in 2D. Thus, the 3D 
and 4D results support the massive solution of the SDEs 
[3, 4], while the 2D case has a scaling behavior Then, a 
natural question is why the 2D case is different from the 
3D and 4D ones. At the moment, a possible answer to 
this question has only been presented in [5]. 

Recently, three works [6, 7, 8] have allegedly shown 
evidence of the scaling IR behavior also in 3D and 4D. 
Here, we will comment on these three works. 
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pact) definition and the (non-compact) stereographic pro- 
jection [10]. Their main conclusions are: "We further- 
more demonstrate that the massive branch observed for 
a^q^ < 1 does depend on the lattice definition of the 
gluon fields, and that it is thus not unambiguously de- 
fined.... One might still hope that this ambiguity will go 
away at non-zero j3 in the scaling limit. While this is 
true at large momenta, we demonstrate.. .that the ambigu- 
ity is still present in the low-momentum region, at least 
for commonly used values of the lattice coupling such 
as j3 = 2.3 or P — 2.5 in SU {2). ...The scaling proper- 
ties such as exponent and coupling, on the other hand, 
appear to be robust under variations of the discretization 
of the gauge fields. ..This emphasizes the importance of 
understanding any discretization ambiguity of the asso- 
ciated gluon mass, before concluding that this mass is 
now firmly established." However, nowhere in Ref. [6] 
are data at j3 =2.5 shown. On the other hand, data for 
a lattice volume 32"* at j3 = 2.5 in the SU(2) case are 
presented in Ref. [10] for the two propagators, using the 
standard discretization and the stereographic projection. 
The conclusion of [10] is that "...there are hardly any dif- 
ferences between the propagators obtained in each case". 
Thus, referring to the last sentence reported above from 
Ref. [6], there are no discretization ambiguities in the 
evaluation of these propagators and the existence of a 
gluon mass is now firmly established. 



We have already criticized Ref. [6] in our work 
[9]. Since that criticism has not been answered, we 
will repeat it here. The authors of Ref. [6] study the 
Landau-gauge gluon and ghost propagators in the strong- 
coupling limit of pure SU(2) lattice gauge theory. These 
propagators are evaluated using different discretizations 
of the gluon field and, in particular, the standard (com- 



The Absolute Landau Gauge 

Ref. [7] considers the absolute Landau gauge, i.e. con- 
figurations belonging to the fundamental modular region 
A. This approach, however, cannot yield an IR-enhanced 
ghost propagator in 3D or in 4D. Actually, restricting 



the configuration space to the region A makes the ghost 
propagator even less singular [11]. This can be seen, in- 
deed, also in Figures 5 and 12 of Ref. [7]. The author 
tries to explain these results, which clearly go against the 
scaling solution, by saying that "The reason for this be- 
havior of the ghost propagator. .may be connected to the 
volume evolution of the first Gribov region and the fun- 
damental modular region. ...The combined effect of the 
precise shape of the low-eigenvalue spectrum and a di- 
verging normalization of the eigenstates could be suffi- 
cient to provide a more infrared divergent ghost propaga- 
tor in the infinite-volume and continuum limits in abso- 
lute Landau gauge than in minimal Landau gauge." Thus, 
simulations in the absolute Landau gauge should agree 
with the scaling solution in the infinite-volume limit (for 
D =3,4) due to a hypothetical diverging contribution 
of the eigenstates of the Faddeev-Popov operator (Note 
that a possible way of quantifying this sentence would 
be to prove that the lower bound of the ghost propaga- 
tor, introduced in Ref. [12], blows up sufficiently fast 
in the infinite-volume limit.) Moreover, that this effect 
should be important in the absolute Landau gauge and 
not in the minimal Landau gauge remains a mystery to 
us, considering that any configuration belonging to the 
absolute Landau gauge is also a configuration of the min- 
imal Landau gauge. The author also adds "A final proof 
is, of course, only that in the absolute Landau gauge at 
sufficiently large volume the ghost propagator would be 
more singular than in the minimal Landau gauge. The 
volume dependence of the propagator in both gauges 
found here is as expected if this is the case. ...Hence, at 
the current point it seems more appropriate to compare 
the lattice results in absolute Landau gauge, rather than 
in minimal Landau gauge, with those from functional re- 
sults which exhibit a scaling behavior in the far infrared." 
These statements — which are somewhat sibylline, since 
the data do not show an IR-enhanced ghost propagator 
— may be the reason why Ref. [7] is (wrongly) cited 
as a numerical verification of the scaling behavior for D 
>2. 



The 2D Case 

Let us note that if the massive behavior observed in 
3D and in 4D could be related to discretization effects, 
as suggested by Ref. [6], or to Gribov-copy effects, as 
reported in [7], then these effects should also be present 
for D = 2 and one should not find in this case a scaling 
behavior In this respect, still in Ref. [7], the author 
makes the following prediction: "A consequence of this 
scenario is that it should be expected that also in two 
dimensions, for sufficiently large volumes and number 
of Gribov copies, an infrared finite gluon propagator is 
obtained in the minimal Landau gauge." We checked this 
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FIGURE 1. Plot of D{0)/V, together with its upper and 
lower bounds [13], as a function of the inverse lattice side 1 /N. 
In the three cases we get a behavior l/N'' with e = 2.67(2). 



prediction by evaluating the gluon propagator for lattice 
volumes up to 2560^ at j3 = 10 (i.e. with a lattice size 
L Ri 460 /ot). In Fig. 1 we plot the volume dependence 
of D{0) /V, together with the upper and lower bounds 
introduced in Ref. [13]. The three sets of data clearly 
extrapolate to zero faster than 1 /V, implying D{Q) ~ 
in 2D. 



The B Gauges 

After considering the absolute Landau gauge in Ref. 
[7], a new set of gauges — called B gauges — was in- 
troduced by the same author [8]. In this case one looks 
along each orbit for a transverse configuration that yields 
a given value B for the ghost dressing function Dq{p) = 
p^G{p) at the smallest non-zero momentum /?„„■„. This 
definition does not solve the Gribov ambiguity [8]. More- 
over, in order to find an IR-enhanced ghost propagator 
one needs to favor configurations closer to the first Gri- 
bov horizon dH. This is the opposite of what is done in 
the absolute Landau gauge, where one favors configura- 
tions well inside the first Gribov region £2 [14]. Thus, 
if the B gauges should produce the scaling solution, "it 
could well be that. . .the absolute Landau gauge is not con- 
nected to a scaling behavior" [8], in disagreement with 
Ref. [7]. 

The main result of this approach is that the ghost prop- 
agator is strongly affected by the choice of configura- 
tion on each orbit, in such a way that its values are en- 
closed in a "corridor". In particular, in 3D and 4D the 
upper bound of this corridor "is strongly increasing with 
volume". At the same time, the gluon propagator seems 



to be B-independent and we should have D(0) > in 
the infinite-volume limit. Thus, the only scaling solution 
that can be obtained with the B gauges seems to be the 
one corresponding to a critical exponent k: = 1/2, which 
was never the preferred value in scaling-solution works. 
Moreover, if the infrared exponent is 1/2 then in 4D one 
should have Dg(p) ~ 1 Iv- Since 1 / pmm ~ i-, the upper 
bound of the corridor should grow at least as fast as the 
lattice size L, in order to support the scaling solution. 
One can verify that this is not the case with the 4D data 
presented in Fig. 3 of [8] (the curve should be hyperbolic 
as a function of 1 /L). 

Let us note that one of the motivations for the intro- 
duction of the B gauges is the possible relation with the 
one-parameter family of solutions obtained by functional 
methods [3, 15|. In this respect one should stress, how- 
ever, that the B gauges are related to different Gribov 
copies on each orbit. On the other hand, the configura- 
tion space is not encoded in the SDEs and this informa- 
tion has to be put in by hand. This can be done in simple 
cases [16], if all Gribov copies are known, but nobody 
knows how to do it in a realistic case. Thus, this relation 
seems at the moment quite accidental. Even more fanci- 
ful seems to us the hypothetical connection between the 
Kugo-Ojima (KO) approach [1] and a (possible) scahng 
solution obtained using B gauges. Indeed, this connec- 
tion requires "subtle cancellation" [8], since one has to 
relate an average over all Gribov copies to results ob- 
tained by selecting specific copies inside the first Gribov 
region Q. In our opinion, the lack of BRST invariance 
when the functional space is restricted to H [17] obscures 
the relation of the GZ approach with the KO criterion and 
the analogies between these two approaches seem to be, 
at the moment, a questionable coincidence. 

Finally, several questions should be answered before 
discussing in detail the results obtained using B gauges. 
For example, it is well known that some Gribov copies on 
the lattice are just lattice artifacts. Thus, by using the B 
gauges, aren't we just probing these artifacts? This may 
explain the over-scaling observed in [18]. Also, it seems 
very difficult to control the infinite-volume limit of the 
corridor and, as pointed in [8], "it cannot be excluded 
that the corridor closes again at much larger volumes". 
This seems indeed possible since, for very large lattice 
volumes, all the orbits should come very close to the 
boundary dQ. and one can expect smaller Gribov-copy 
effects [19]. 

CONCLUSIONS: THE GHOST FACTORY 

We believe that, in order to understand the results ob- 
tained in minimal Landau gauge using numerical simu- 
lations, the first question to be answered is: why is the 
2D case different? One could also ask: can we test nu- 



merically the explanation presented in |5|? A clear an- 
swer to these questions probably requires new ideas and 
better data (especially in the ghost sector). Unproven hy- 
potheses and happy coincidences should on the contrary 
be treated with great caution. 

We recently installed at IFSC-USP a new machine 
with 18 CPUs Intel quadcore Xeon 2.40GHz (with In- 
finiBand network and a total of 216 GB of memory) 
and 8 NVIDIA Tesla SI 070 boards (500 Series), each 
with 960 cores and 16 GB of memory. The peak per- 
formance of the 8 Tesla boards is estimated in about 
2.8 Tflops in double precision and 33 Tflops in single 
precision. This machine will be used mainly for studies 
of Green's functions in different gauges (Landau, Feyn- 
man and Coulomb) for various SU(A^t) gauge groups. In 
particular, the GPUs will be used for the inversion of 
the Faddeev-Popov matrix using conjugate gradient. This 
computer will allow us to perform an extensive study of 
the ghost sector. We believe that this new ammo will help 
us clarify the issues addressed above. 
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